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We study a ferromagnetic tendency in the two-dimensional Hubbard model near van Hove 
filling by using a functional renormalization-group method. We compute temperature depen- 
dences of magnetic susceptibilities including incommensurate magnetism. The ferromagnetic 
tendency is found to occur in a dome-shaped region around van Hove filling with an asym- 
metric property: incommensurate magnetism is favored near the edge of the dome above 
van Hove filling whereas a first-order-like transition to the ferromagnetic ground state is ex- 
pected below van Hove filling. The dome-shaped phase diagram is well captured in the Stoner 
theory by invoking a smaller Coulomb interaction. Triplet p-wave superconductivity tends 
to develop at low temperatures inside the dome and extends more than the ferromagnetic 
region above van Hove filling. 
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The stability of itinerant-electron ferromagnetism is an important problem of correlated 
d-electron systems.^' Stoner proposed a simple approach to treat the ferromagnetic insta- 
bility, introducing an averaged magnetic field of electrons at each lattice site. His theory is, 
however, insufficient to describe adequately even magnets with small magnetic moments, es- 
pecially at finite temperatures. The transition temperatures appear too high in comparison 
with experimental data and the temperature dependence of the magnetic susceptibility is not 
described correctly. 

Dzyaloshinskii, Kondratenko^^ and Moriya"^^ proposed the spin-fluctuation theory consid- 
ering a contribution of paramagnons to thermodynamic properties. They improved the Stoner 
theory, especially for finite temperature properties of weak and nearly ferromagnetic systems. 
The condition for ferromagnetic order in the ground state is also corrected by their theory. 
In the regime of strong Coulomb repulsion the stability of ferromagnetism was considered in 
the pioneering studies by Nagaoka^^ and Roth,^) which showed the existence of a saturated 
ferromagnetic region close to half filling and non-saturated ferromagnetism further away from 
it. The region of the saturated ferromagnetism was investigated in more details within various 
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approximations in numerous works (see discussion and references in Ref7^). 

Two-dimensional systems offer an interesting aspect to study the stability condition of a 
ferromagnetic ground state and its finite-temperature properties. Because of van Hove singu- 
larities in the density of states ferromagnetism can be realized within a one-band model even 
for moderately small Coulomb repulsion U.^^ This allows to apply perturbative techniques 
to this problem to verify and understand results of other approaches. At finite temperatures, 
the continuous symmetry breaking is absent in two dimensions due to fluctuation effects, re- 
sulting in fulfillment of the Mermin- Wagner theorem; the Curie- Weiss law is expected to be 
accordingly modified. 

Another peculiarity of two-dimensional systems is that the ferromagnetic order in two 
dimensions can be destabilized not only by quantum magnetic fluctuations, e.g., via spin- 
wave excitations, which are nearly commensurate, but also by incommensurate fluctuations. 
The importance of incommensurate magnetic fluctuations at weak and moderate Coulomb 
interaction was recently emphasized within the quasistatic approach,^' mean-field theory,-*^^) 
and a renormalization-group approach. These approaches showed that in a large part of the 
phase diagram spanned by the electron density n and Coulomb interaction U the ferromagnetic 
order is replaced by an incommensurate one. 

Therefore, it is important to consider systematically the effect of commensurate and incom- 
mensurate magnetic fluctuations on the ground-state phase diagram and finite temperature 
properties of ferromagnetism in two dimensions. In an attempt to investigate this problem, 
we use a functional renormalization-group (fRG) approach in the symmetric phase^^^^^^ and 
study temperature dependences of magnetic susceptibilities. By extrapolating these results 
to the limit of zero temperature, one can obtain expected properties such as possible insta- 
bilities of the ground state. This procedure was recently applied to studying a possibility of 
antiferromagnetic and superconducting instabilities near half filling. 

Model and method. We consider the two-dimensional (2D) t-t' Hubbard model = 
H-{n- 4t')N with 

H = -Y1 tijcj^Cja + rii^rii^ , (1) 

ija i 

where tij = t for nearest neighbor (nn) sites i, j, and tij = —t' for next-nn sites (t, t' > 0) on 
a square lattice; c|^(cjcr) creates (annihilates) an electron with spin (J dbt i site; tii^j — c|^Qcr 
and = ^ 

i(j f^ia- For convenience we have shifted the chemical potential /i by At' . We choose 
t' /t = 0.45 as a typical value, for which the ferromagnetic instability is favored.^^'^^'^^^ We 
employ the fRG approach for the one-particle irreducible generating functional and choose the 
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temperature T as a natural cutoff parameter as proposed in Ref.^^^ Neglecting the frequency 
dependence of interaction vertices, the RG differential equation for the interaction vertex 
Vt = y(ki, k2, k3, k4) (with the momenta kj , which are supposed to fulfill the conservation 
law ki + k2 = ka + k4) has the form^^"^^) 

dVr dLpp dLph 

= —Vt o — —oVt + Vto — —oVr, (2) 

dT ^ dr ^ ^ dT ^ ' ^ ^ 

where o is a short notation for summations over intermediate momenta and spins, -^ph,pp 

stand for particle-hole and particle-particle bubbles. The ferromagnetic (FM), triplet p-wave 

superconducting (pSC) and incommensurate magnetic (Q) susceptibilities can be calculated 

as 

'^Xm Y^-r> m T?"i dLph,pp(k ; q^) 

-J^ = 2^ '^-k'.ik'+q^'^-k'.ik'+q^ -T^ > 

k' 

d^k^zfck-q^ Sr-^m pT/, WN^^^ph^PpO^lSlrrO 

==F Z^''*^k',±k'+q^iml^>*^J ^ > 

k' 

where the three-point vertices T^J^^' describe the propagation of an electron in a static external 
field, m denotes a type of instability, qpM.pSC = and qq = Q; upper signs and ph correspond 
to the magnetic instabilities, lower signs and pp to the superconducting instability; 



r^(k,k') 



l^(k,k',k' + q^,k-q^) m = FMorQ, 
F(k,-k,k',-k') m = pSC. 



Eqs. (2) and (3) are solved with the initial conditions VTo(ki, k2, ks, k4) = U, 'R-i^k+q ~ 
/k and Xm = 0; the initial temperature is chosen as large as Tq = lO^t. The function /k 
belongs to one of the irreducible representations of the point group of the square lattice, 
i-e., /k = 1 for the magnetic instabilities and /k = sin A;,i.,y/^ for the pSC, with A being a 
normalization coefficient. We also discrctizc the momentum space in Np = 48 patches using 
the same patching scheme as in Ref."*^^-* This reduces the intcgro-differential equations (2) and 
(3) to a set of 5824 differential equations, which were solved mimcrically. In the present paper 
we perform the renormalization-group analysis down to the temperature T^q^, where vertices 
reach a maximal value; we choose Vmax = 18t. 

To characterize the strength of fluctuations of the order parameter m, we introduce the 
temperature which is defined by the condition x^^(X^) = 0. Here x~^(r) is the inverse 
susceptibility of the order parameter m, analytically extrapolated to the region T < T^q^ 
(see Ref.^^) for details). We interpret as a crossover temperature to the regime of strong 
correlations rather than a phase transition temperature, since the finite temperature tran- 
sition is prohibited by the Mermin- Wagner theorem, inducing non-analytical corrections to 
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susceptibilities, e.g., Xfmq(^) e^^^^-'^/-^ at low T as a consequence of an exponentially 
large correlation length (see, e.g., Ref.-*^^)). 

To compare the obtained results with the Stoner theory we also study Hamiltonian 
(1) in the mean-field approximation, by decoupling ni-^rii^ ("t)'^i4- + ^it('^-l-) ~ ('^t)(%)- 
The second-order transition temperature (Curie temperature) is obtained by the condition 
Xp^(r^^) = 0. To search for possible first-order transitions, we also compute the Landau 
free energy as a function of magnetization for fixed chemical potentials. 

Results. In Fig. 1(a) we present the results of the fRG approach for the inverse magnetic 
susceptibilities, obtained at U = At and different fillings above van Hove filling n = 0.465 (/Lt = 
0). Above the crossover temperature the susceptibilities follow approximately the Curie- Weiss 
law with Curie temperatures replaced by Tpj^. With increasing ji, the inverse susceptibilities 
increase, indicating that ferromagnetic fluctuations are weakened. Consequently, the crossover 
temperature Tpj^ decreases and vanishes at ~ 0.052t, implying the quantum phase transition 
from the ferromagnetic to paramagnetic phase. This transition is, however, likely preempted 
by entering into the phase with strong incommensurate fluctuations, since we conflrmed by 
calculations of Xq, which yield results similar to Fig. 1(a), that in the range jj, G (0.05, 0.055)t 
a condition Tq > Tpj^ is fulflUed for some of the vectors Q = (Q, Q) with small Q. Therefore, 
for the ground state above van Hove flUing one can expect two successive quantum phase 
transitions: from the ferromagnetic to an incommensurate and then to the paramagnetic 
phase with increasing ^. A direct phase transition from the ferromagnetic to paramagnetic 
phase having strong incommensurate magnetic fluctuations is yet also possible. Further away 
from van Hove filling we obtain a non-monotonic temperature dependence of the inverse 
susceptibility with a minimum, which is followed by a maximum for /x G (0.065, 0.085)t. 
While the minima of the inverse susceptibility away from van Hove filling are produced by 
thermal excitations of states near van Hove singularity, an interpretation of the maxima 
which possibly yield a 'reappearance' of magnetic order far away from van Hove filling is not 
so straightforward. We return to the discussion of this peculiarity below. Even further away 
from van Hove filling [n > 0.09i) the flow can be continued down to very low temperatures, 
such that T^^ « 0. 

Let us consider now the situation below van Hove filling, which is distinct from the situa- 
tion above van Hove filling. The obtained temperature dependence of the inverse susceptibil- 
ities [see Fig. 1(b)] shows larger bending than that for the Fermi level above van Hove filling; 
the inverse susceptibilities show deviations from the Curie- Weiss law already at T > Tpj^. 
Moreover, the obtained crossover temperature Tpj^ sharply drops to zero at /x ~ — 0.065t. The 
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Fig. 1. (Color online) Temperature dependences of the inverse magnetic susceptibility at U = 4t and 
different values of the chemical potential (marked by numbers) for the Fermi level above (a) and 
below (b) van Hove energy. Dots denote the temperature T^q", at which the fRG flow is stopped; 
note that T^]^^ ~ for /i = O.lOi. Below T^q", the susceptibility is analytically extrapolated. 

precise form of the dependence of Tp^{^) in this region depends on the details of analyti- 
cal extrapolation. One can obtain either discontinuous or sharp continuous transition with 
vanishing T^^. Since Tpj^ should be interpreted as a crossover temperature characterizing 
magnetic properties of the ground state in a 2D system, this necessarily implies a sharp change 
of those around fi ~ — 0.065t, in particular a possibility of a first-order phase transition from 
the ferromagnetic to paramagnetic ground state. Studying the susceptibilities of incommen- 
surate magnetic orders shows that we always obtain Tq < below van Hove filling, i.e., 
incommensurate magnetic fluctuations are not expected to change our conclusion. 

The resulting finite-temperature phase diagram in the T-/i variables is shown in Fig. 2. 
Tpjyj forms a dome-shaped line, which is asymmetric with respect to the van Hove energy 
(/i = 0). The incommensurate magnetic tendency occurs near the edge of the dome above the 
van Hove energy whereas such a feature is not seen below the van Hove energy and instead 
a first-order-like transition to (commensurate) ferromagnetism is expected as a ground state. 
In Fig. 2, we also plot T^q^{^), which is usually interpreted as a crossover temperature to an 
ordering tendency, ^^"^^^ by rescaling both temperature and /x by a factor of 1.6. The line of 
almost coincides with that of T^^, which suggests that our extrapolation procedure to 
obtain Tpj^ is performed in a reasonable way. Employing the same extrapolation procedure 
for the non-monotonic dependence of the inverse susceptibilities seen in Fig. 1(a) around 
/X ~ 0.07t, we obtain the 'second', tiny, ferromagnetic region as shown in Fig. 2. We are 
however not aware of the physical explanation of this possible 'reappearance' of ferromagnetic 
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order, and leave it for future studies. 

Since triplet superconductivity is expected around a ferromagnetic state/^^ we also com- 
pute XpSC as a function of T in the same fashion of Fig. 1. In Fig. 2 we plot T*gQ, which, similar 
to previous results, is smaller than Tpj^ around van Hove filling and survives even far above 
van Hove filling, where p-wave superconductivity is therefore expected in the ground state. 
On the other hand, we do not find an appreciable T*gQ away from the ferromagnetic phase 
below van Hove filling. In the region /i G (—0.07, 0.05)t the coexistence of ferromagnetism and 
p-wave superconductivity is possible, which is a subject of future studies. 

0.06 



0.04 
0.02 


-0.15 -0.1 -0.05 0.05 0.1 0.15 

Fig. 2. (Color online) Phase diagram obtained in the fRG study for U = 4Hn the plane of /i and T. 
PM denotes the paramagnetic phase. 

We now compare the fRG results to those obtained in the Stoner theory, which considers 
only a (commensurate) ferromagnetic instability. Since the Stoner theory predicts much higher 
transition temperatures and a much broader concentration range of ferromagnetism, we 
consider smaller U = U"^^, which is chosen such that max^ {fi;U^^) = max^ Tpj^(^; [/). 
For U = At we obtain U^^ ~ l.7t. This renormalization of U is due to fluctuations, mainly 
caused by particle-particle scattering processes, similar to the original idea by Kanamori.^*^) In 
Fig. 3, we see that not only the height, but also the position and the width of the ferromagnetic 
region are in good agreement with those in the fRG approach. Furthermore, below van Hove 
filling, the possibility of the first-order-phase transition as a function of fi (Fig. 2) within 
the fRG is also in agreement with the results of 'renormalized' Stoner theory, since such a 
transition is transformed to a phase separation in terms of electronic density. ^"^'^^^ Above 
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Fig. 3. (Color online) Comparison of the fRG phase diagram in the n-T coordinates with the Stoner 
theory and with the ground-state results of mean-field approach including incommensurate mag- 
netic orders^^-* and T-matrix approximation.^^) PS denotes a region of phase separation. 



van Hove filling the renormalized Stoner theory yields a region of phase separation, while 
the fRG predicts the strong tendency of incommensurate magnetism in the corresponding 
region. Except for this difference, it is remarkable that the ferromagnetic tendency obtained 
in the fRG is captured very well in the renormalized Stoner theory with U^^, which is n and 
T-independent . 

Considering incommensurate magnetic order in the mean-field theory yields a variety of 
the ground states in the plane of n and U , as revealed in Ref.^^^ In particular, the mean-field 
theory^^^ already predicts incommensurate magnetism above van Hove filling and (commen- 
surate) ferromagnetism roughly below it, for U = At the ground state changes at n = 0.55 
through a second-order transition. Naturally, the mean-field treatment of Ref.^^-* predicts mag- 
netic order in a very broad region (0 < n < 0.9), which is obviously an overestimation, typical 
for a mean-field theory. In particular, a broad incommensurate region (0.55 < n < 0.9) is 
replaced by a very narrow region with strong incommensurate fluctuations in the fRG ap- 
proach. Finally we find that the positions of the quantum phase transitions obtained in the 
fRG approach show reasonable agreement with the results of T-matrix approach. The T- 
matrix approach, however, misses the importance of incommensurate magnetic fluctuations 
above van Hove filling. 

Conclusion. We have studied finite temperature phase diagrams and possible ground-state 
properties of the one-band Hubbard model at different chemical potentials or fillings. The re- 
sults are compared to those in the Stoner theory, including a possibility of incommensurate 
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magnetic order. We have found that a ferromagnetic tendency occurs in a dome-shaped region 
around van Hove filhng in an asymmetric way: incommensurate magnetism is favored near 
the edge of the dome above van Hove fiUing, whereas a first-order-hke transition to a (com- 
mensurate) ferromagnetic ground state is expected below van Hove filling. The verification of 
the first-order transition however requires further development of the present fRG approach, 
e.g., its extension to the symmetry-broken phase. At finite temperatures we observe devia- 
tions from the Curie- Weiss behaviour of susceptibilities, which are mostly pronounced below 
van Hove filling and seen already above the temperature of the crossover to the regime of 
strong magnetic correlations. In agreement with previous studies, a ferromagnetic tendency 
is accompanied by development of triplet superconductivity at low temperatures. While our 
study indicates a pure p-wave superconducting state far above van Hove filling, we cannot 
address whether the ground state is the coexistence of ferromagnetism and triplet supercon- 
ductivity in a region where a ferromagnetic tendency also occurs. Given that the coexistence is 
actually observed for various U-based compounds such as 1)062,^^-* URhGe^^^ and UCoGe,^^^ 
it is an interesting subject to elucidate a possible coexistence of ferromagnetism and ^?-wave 
superconductivity in the ground state of the Hubbard model near van Hove filling. 
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